Strong-coupling study of the Gribov ambiguity in lattice Landau gauge 



O 

Oh. 
< 

m 



Oh' 



(N 
> 
m 
o 

(N 

ON 

O 



X 



Axel Maas,^ Jan M. Pawlowski,^''^ Daniel Spielmann,^' Andre Sternbeck,*' ^ and Lorenz von Smekal® 

'institut fiir Physik, Karl-Franzens Universitdt Graz, Universitdtsplatz 5, A-8010 Graz, Austria 
^Institut fiir Theoretische Physik, Universitdt Heidelberg, Philosophenweg 16, D-69120 Heidelberg, Germany 
ExtreMe Matter Institute EMMI, GSI Helmholtzzentrum fiir 
Schwerionenforschung, Planckstr. 1, D-64291 Darmstadt, Germany 
'^Centre for the Subatomic Structure of Matter, University of Adelaide, SA 5005, Australia 
''Institut fiir Theoretische Physik, Universitdt Regensburg, D-93040 Regensburg, Germany 
''Institut fiir Kernphysik, Technische Universitdt Darmstadt, Schlossgartenstr. 9, D-64289 Darmstadt, Germany 

We study the strong-coupling limit /3 = of lattice SU(2) Landau gauge Yang-Mills theory. In 
this limit the lattice spacing is infinite, and thus all momenta in physical units are infinitesimally 
small. Hence, the infrared behavior can be assessed at sufflciently large lattice momenta. Our results 
show that at the lattice volumes used here, the Gribov ambiguity has an enormous effect on the 
ghost propagator in all dimensions. This underlines the severity of the Gribov problem and calls for 
refined studies also at finite /3. In turn, the gluon propagator only mildly depends on the Gribov 
ambiguity. 
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I. INTRODUCTION 

The study of the long-range behavior of QCD Green's 
functions is one path towards the understanding of con- 
finement [1;, 2]. In particular, the confining properties of 
QCD are directly reflected in the behavior of the Landau 
gauge gluon and ghost propagators at asymptotically low 
momenta [3j. Furthermore, confinement, as signaled by 
an unbroken center symmetry, can also be shown to be a 
direct consequence of these properties Q ■ 

However, implementing Landau gauge as in the pertur- 
bative regime is not sufficient to uniquely fix the infrared 
behavior of Green's functions because of the Gribov- 
Singer ambiguity [F, Using functional methods, one 
finds a one-parameter family of solutions which can be 
classified by the zero momentum value of the ghost prop- 
agator 0, |3| ■ These can be obtained by specifying a non- 
perturbative renormalization condition [7|. 

The one-parameter family of solutions contains the 
scaling solution as a limiting special case 0, which is 
characterized by an infrared divergent ghost dressing 
function, and an infrared suppressed and likely vanish- 
ing gluon propagator, see e.g. [Hlol-fTTj. In contrast, the 
other members of the family with an infrared finite ghost 
dressing function necessarily come along with an infrared 
finite gluon propagator, see e.g. [tI. IsI. fl8l - l26j . These are 
called the decoupling solutions Both terms are more 
precisely defined in Section HIl 

These Green's functions have also been determined on 
the lattice. By now it is well established that standard 
lattice gauge fixings lead to a decoupling- type behavior, 
see e.g. ^27- 37], though not unanimously [H, HI]. More- 
over, investi gat ions of the Gribov-Singer ambiguity in 
this setting, Issl-iaTlliol-lisj . so far only led to rel- 

atively small, quantitative effects. 

However, recently alternative gauge fixing procedures 
have been devised that directly make use of the Gribov- 



Singer ambiguity [4y|. Here, one selects different non- 
perturbative completions of the Landau gauge by con- 
stricting the ghost propagator in the zero momentum 
limit to be as close as possible to a fixed value. Such a 
procedure samples Gribov copies differently than a stan- 
dard gauge fixing. It introduces an additional parame- 
ter which plays a role analogous to the non-perturbative 
renormalization condition of the continuum studies. The 
interpretation here is that this free parameter is there- 
fore linked to the residual gauge freedom of the Landau 
gauge. 

Note that the qualitatively different infrared solutions, 
decoupling and scaling, signal different realizations of 
the global part of the gauge fixing. At the level of 
gauge- invariant observables they have to be equivalent. 
One particular difference is the structure of the Becchi- 
Rouet-Stora-Tyutin (BRST) symmetry. If a given non- 
perturbative completion of the Landau gauge is to have 
a global BRST symmetry then the Kugo-Ojima confine- 
ment criterion requires that the corresponding Landau 
gauge ghost propagator is infrared enhanced, correspond- 
ing to the scaling solution 0, E^]. On the other hand, 
in the decoupling case, the infrared ghost propagator is 
essentially massless. In standard Faddeev-Popov the- 
ory, the corresponding massless asymptotic ghost states, 
albeit being unphysical, would then contribute to the 
global gauge charges and lead to a spontaneous breaking 
of the global gauge symmetry much like what one expects 
for theories with a Higgs mechanism. When there is no 
Higgs mechanism, and thus no physical vector field in the 
gauge boson sector, this breaking is most likely artificial 
and due to a non-perturbative breakdown of BRST sym- 
metry itself (23l. [24 lisl - lssj l . At present, the only solution 
which is consistent with an unbroken BRST symmetry 
is thus the scaling solution. Note that other BRST for- 
mulations might exist, in which the massless ghosts of 
the decoupling solutions would not automatically lead 
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to BRST breaking. However, such a formulation is still 
lacking, but investigations along the lines of [sil appear 
promising. 

On the lattice, the only way that is currently known to 
define BRST symmetry in presence of Gribov copies and 
thus fully non-perturbatively is based on defining gauge 
fields by stereographic projection. This avoids the per- 
fect cancelation of Gribov copies and the Neuberger 0/0 
problem of lattice BRST, while it nevertheless reduces 
to the standard BRST symmetry in the continuum limit, 

MM- 

As noted above, for studying the long-range behavior 
of the Landau gauge gluon and ghost propagators on the 
lattice, it is useful to employ the strong-coupling limit 
(i^ - H . HH, 113 • In our present work we use the strong- 
coupling limit /? — > for pure SU(2) lattice gauge theory 
in d = 2 and 3 dimensions. This limit can be interpreted 
as the limit of infinite lattice spacing, a — > oo, at a fixed 
physical scale as set, e.g., by the string tension. Alter- 
natively, when considering all momenta in lattice units, 
q ~ 1/a, the same limit can also be interpreted as the 
hypothetical limit in which the physical scale is sent to 
infinity, i.e., formally as Aqcd ^ oo. Both interpreta- 
tions are equivalent. They are based on taking the limit 
of the function a(/3) for /3 — 0. Other ways to assign 
a scale to the theory at (3 — are possible [6ll, [gl] , but 
they require different prescriptions. For the limit /3 — >■ 
adopted here, all momenta and masses in lattice units 
1/a are infinitely small relative to the physical scale of 
the theory which is precisely what we need for an analy- 
sis of the asymptotic infrared behavior of its correlation 
functions. It is therefore best suited to isolate this be- 
havior from finite- volume effects. 

In recent studies (il - lil. l6lj | the /3 -)■ limit was first 
investigated only in four dimensions |42h44| and then to 
some extent also in lower dimensions [6l|. The results 
of both studies were compatible with a scaling behav- 
ior for the gluon propagator for intermediate to large 
momenta and a decoupling behavior at small momenta, 
where 'large' and 'small' in the strong-coupling limit al- 
ways refers to lattice units. It has moreover been shown 
in [42I - I43 I that the mass parameter related to the low- 
momentum decoupling tail can be changed by choosing 
different definitions of the gluon field, which should be 
equivalent in the continuum limit. For non- vanishing lat- 
tice spacing a, however, these correspond to different lat- 
tice implementations of the Landau gauge [5&] . If some 
of these differences remain in the continuum limit that 
would strengthen the hypothesis 0, ^ residual 

gauge freedom related to the one-parameter family of so- 
lutions observed in the continuum. In any case, the lat- 
tice discretization dependence adds to the puzzle of the 
global properties of Landau gauge Yang-Mills theory. 

The interpretation of the data for the ghost propaga- 
tor has always been even less conclusive, and did not 
allow a direct determination of local scaling coefficients. 
This led to a more indirect interpretation of the data. 
In [4344^ it was concluded that the ghost data at larger 



lattice momenta could be consistent with the scaling be- 
havior observed in the gluon propagator. However, it 
was not possible to extract a firm value for the scaling 
exponent from the ghost data for any reasonably wide 
range of lattice momenta. In turn, the authors of [6"T| 
assessed the evidence of a logarithmic momentum de- 
pendence and argued in favor of decoupling at all lattice 
momenta for Landau gauge Yang-Mills theory in three 
and four dimensions in the strong-coupling limit. It was 
also frequently claimed that the scaling solution is real- 
ized in two dimensions in agreement with the arguments 
in [eil and the lattice data of [H, 113, Ull . 

In this paper we present results with much improved 
statistical accuracy in two and three dimensions. We also 
provide a more detailed study of the systematic uncer- 
tainties, see Appendices [K\ and [Bl and of the relevance of 
the Gribov ambiguity, see Section HIIl and parts of Sec.lIVI 

Our results are presented in Sec. |TVl For th e g luon 
propagator we confirm the scaling results of [i^ - H l6l| 
at intermediate lattice momenta. In particular, we obtain 
fairly stable scaling exponents also with the local analy- 
sis of [6l| . The values for the exponents tend to decrease 
somewhat at large lattice momenta, and as in the previ- 
ous studies we observe decoupling at small momenta, at 
least in three dimensions. 

A qualitatively new and perhaps surprising behavior is 
seen in the ghost propagator which shows a very strong 
sensitivity to Gribov copies at low momenta. This en- 
tails that we still have to work on a better understand- 
ing of the global properties of Landau gauge QCD. In 
particular our analysis shows that a proper treatment of 
Gribov copies in the thermodynamical limit of Landau 
gauge lattice Yang-Mills theory has not been achieved 
yet. Subject to the sampling of Gribov copies the ghost 
propagator can be tuned to also show a scaling behavior 
for intermediate momenta. The general trend is that even 
a drastic over-scaling is obtained if enough gauge copies 
are taken into account in the search for the maximally 
achievable infrared enhancement. In fact, the possibility 
of such an over-scaling has also been observed in con- 
tinuum Landau gauge Yang-Mills theory on the torus 
[5^ . Here, we moreover demonstrate that this peculiar 
behavior arises already in two dimensions where infrared 
scaling was supposed to be quite well established. We 
find that the two-dimensional ghost propagator shows 
the same apparently unconstrained infrared behavior as 
that in higher dimensions, and that there is thus no qual- 
itative difference between two and higher dimensions as 
long as sufficiently many Gribov copies are present, as 
has been suspected [45, 46]. 

In the summary in Section |V] we are led to conclude 
that our analysis clearly shows the direct relation be- 
tween the ambiguities in the ghost propagator at low and 
intermediate momenta and the Gribov ambiguity. In our 
opinion this ambiguity is much more than an artificial 
problem of this unphysical limit. In fact, Gribov copies 
at finite lattice coupling l3 will remain to affect the in- 
frared behavior of Green's functions as long as they are 
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present in the strong-coupling limit. To resolve it there 
will certainly be more cost-efficient than to try to reach 
the continuum limit in a sufficiently large physical vol- 
ume by brute force, though this is attempted as well [67j . 



II. INFRARED SCALING AND DECOUPLING 

The Landau gauge gluon propagator, in (Euclidean) 
momentum space, is parameterized by a single dressing 
function Z, 

^W-^-^k,-^")^, (1) 



p. 



p. 



and the ghost propagator by a corresponding dressing 
function G, 



(2) 



For their infrared behavior, i.e., that of Z{p^) and G{p^) 
for — 0, we consider the two possibilities of scaling and 
decoupling. Occasionally, we also use Dgi = Z{p'^)/p'^ 
and Dgh = G{pf/p\ 



A. Scaling 

The prediction of [osj amounts to infrared 

asymptotic forms 



Z{P') 



{p'/K 



QCB) 



, o I 4— d 

\2k.z + ^— 



(3a) 
(3b) 



for — >■ 0, which are both determined by an unique 
critical infrared exponent 



Kz = K , 



(4) 



with (d — 2)/4 < K < d/A. Furthermore it has been 
proven in [3| that the scaling solution is confining for 



K > 



d-3 



(5) 



Under a mild regularity assumption on the ghost-gluon 
vertex [l2l | , the value of this exponent is furthermore ob- 
tained in four dimensions as [ill, 113 



(93 - Vl201)/98 w 0.595 , 



(6) 



which is confinement according to ([5]). The correspond- 
ing values in lower dimensions are n — 1/b {d ^ 2) 
and K ~ 0.3976 (d = 3) which also satisfy The 
conformal nature of this infrared behavior in the pure 
Yang-Mills sector of Landau gauge QCD is evident in 
the generalization to arbitrary gluonic correlations [l^ : 



an uniform infrared limit of one-particle irreducible ver- 
tex functions F™'" with m external gluon legs and n pairs 
of ghost/anti-ghost legs, in four dimensions of the form 



(pVa^ 



qcd; 



\ {n — 7n)K, 



(7) 



when all cx p'^ — ^ 0, i = 1, . . . , 2n -f m. In particular, 
the ghost-gluon vertex is then infrared finite (with n = 
m = 1), in agreement with its STI [69], and the non- 
perturbative running coupling introduced in 0, [l3| via 
the definition 



asip') = '-Z{p^)G\p^) 



(8) 



approaches an infrared fixed point, Ug — )■ etc for p^ ^ in 
d = 4. For arbitrary dimensions, trivial kinematic factors 
lead to an infrared fixed point for the effective coupling 



(9) 



with Q;e(0) = ttc- If the ghost-gluon vertex is regular at 
p^ — 0, its value is in four dimensions [l2| 



Stt r2(K- i)r(4- 2k) 

TVc r2(-K)F(2K- 1) 



— X 0.99 . (10) 



Comparing the infrared scaling behavior of Dyson- 
Schwinger equation (DSE) and functional renormaliza- 
tion group equation (FRGE) for solutions of the form of 
Eqs. it has been shown that in presence of a single 
scale, the QCD scale Aqcd, the solution with the in- 
frared behavior ([4]) and (O, with positive k, is unique 
[l6l . [toI . Iti} . It is nowadays being called the scaling 
solution. 



B. Decoupling 

This uniqueness proof does not rule out, however, the 
possibility of a solution with an infrared finite gluon prop- 
agator, as arising from a transverse gluon mass M , which 
then leads to an essentially free ghost propagator, with 
the free massless-particle singularity at p^ — 0, i.e.. 



Z{p^) ~ //M^ , and G{p^) - const. 



(11) 



for p2 0, H, [T8l-[26j. The constant contribution to 
the zero-momentum gluon propagator, £'(0) = 1/M^, 
thereby necessarily leads to an infrared constant ghost 
renormalization function G. This solution corresponds to 
= 1/2 and k = 0, and it has been proven in [4] that it 
is confining. It does not satisfy the scaling relations ^ 
or ([7]). This is since the transverse gluons decouple for 
momenta p^ <C M^, below the independent second scale 
given by their mass M . It is thus not within the class 
of scaling solutions considered above, and it is termed 
the decoupling solution in contradistinction. The renor- 
malization group invariant ([8|) shows this decoupling as it 
tends to zero in the infrared. Alternatively, we can define 
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a running coupling analogously to that in massive theo- 
ries. By fixing one overall free parameter this definition 
leads to a coupling which even quantitatively resembles 
the fixed-point coupling in the scaling case. This defini- 
tion has been given explicitly in ^7|] and similarly been 
used since then [2^ [2^ . 



III. GAUGE FIXING AND THE GHOST 
PROPAGATOR 

The functional equations of continuum quantum field 
theory, such as DSEs and FRGEs, admit a one-parameter 
family of solutions ^] . The free parameter can be chosen 
as the value of the ghost dressing function at vanishing 
momentum. It fixes an ambiguity in the equations re- 
lated to non-perturbative renormalization and the struc- 
ture of the massless unphysical states in the theory which 
is a priori unknown. The realization of the global gauge 
symmetries and the Kugo-Ojima confinement criterion 
crucially depend on this structure. For that reason func- 
tional equations have been supplemented by the appro- 
priate boundary condition in the earlier literature on the 
scaling solution il2||. The restriction of the gauge field 
configuration space to the first Gribov region provides 
a simple example of boundary conditions which are not 
reflected in the form of the functional equations of the 
theory but need to be imposed in addition . Another 
example is the Kugo-Ojima criterion which cannot be de- 
rived from the functional equations but needs to be added 
as a boundary condition. This then singles out the scal- 
ing solution with an infrared enhancement of the ghost 
propagator. Without this boundary condition, the func- 
tional equations admit both, the scaling solution and the 
decoupling solutions. Since these boundary conditions 
are imposed on unphysical degrees of freedom, they may 
be inconsequential for physical observables. 

Indeed, the Gribov ambiguity and the associated free- 
dom in globally completing the Landau gauge are well 
known to affect particularly the ghost propagator on the 
lattice [13431 IS^M For the volumes accessi- 

ble in current simulations, it appears that this ambiguity 
is indeed closely related to the one-parameter freedom 
in the family of solutions to the functional equations: a 
direct correspondence between the Gribov ambiguity on 
the lattice and the one-parameter family of solutions to 
the functional equations of continuum quantum field the- 
ory is provided by the so-called Landau-B gauges (46l |. 
These gauges correspond to different resolutions of the 
Gribov-Singer ambiguity. They are implemented by first 
selecting a target value, called the B parameter, for the 
ratio of ghost dressing functions evaluated at two widely 
separated momentum scales. An appropriate choice for 
the two momentum scales is the lowest accessible mo- 
mentum on the lattice and some conveniently large mo- 
mentum scale, such as the renormalization scale. This 
choice is made to resemble the non-perturbative freedom 
in choosing a boundary condition for the ghost propaga- 



tor in the continuum studies. Then, on a configuration 
by configuration basis, many Gribov copies are generated 
to select the one where this ratio is closest to the target 
value B. The propagators, or other Green's functions, 
are calculated on these copies from the Monte-Carlo his- 
tory of configurations. 

Any finite number of Gribov copies leads to a finite set 
of proposed B values. In turn, for larger and larger lat- 
tices one can get close to any target value within a certain 
range depending on the lattice volume. We conclude that 
within this range the infrared value of the ghost dressing 
function can be chosen freely. An extreme case is the 
metx-B gauge (46j in which the copy with the maximally 
enhanced ghost propagator is selected. 

In general, it is of course practically impossible to find 
all Gribov copies for a given configuration. This limita- 
tion can cause some residual Gribov noise on top of the 
statistical fiuctuations. In the case of the max-_B gauge 
the result will be a lower bound for the maximally en- 
hanced ghost propagator that is possible in any finite vol- 
ume. In the present work, max-_B gauge is implemented 
similarly to [46| . In order to reduce the effects of the lat- 
tice breaking of rotational invariance, the Euclidean 0{d) 
symmetry, the value of B for a given copy is determined 
by averaging over all representations for the lowest mo- 
mentum, instead of using only one on-axis momentum as 
in [i^. Note also that the standard lattice implemen- 
tation of the Landau gauge, the minimal Landau gauge, 
essentially gives the same results as choosing for B the 
expectation value of the ghost propagator ratio from all 
Gribov copies of all configurations [4g| . 



IV. NUMERICAL RESULTS 

The configurations have been obtained by generating 
randomly j3 — configurations on symmetric lattices of 
size L'^ — (Na)'^. Gauge fixing to the minimal Landau 
gauge in sections [lV Al and llV BI has been perfornied using 
a stochastic overrelaxation method, see e.g. (7^. Fixing 
to the max-B gauge has been discussed in section IIIII 
Propagators are determined using a standard method, 
see, e.g., For the determination of the ghost prop- 

agator we also refer to Appendix A cylinder cut has 
been imposed to select momenta |74| . 

A. Two dimensions 

The gluon propagator for various lattice sizes in d — 2 
is shown in Fig. [1] Its value at vanishing momentum 
continues to decrease with the inverse lattice extension 
(Fig. [2]), and therefore only a lower bound on a possible 
mass at zero momentum can be given. Power-law fits to 
the gluon propagator, according to Eq. ([5al) . yield a value 
of Kz ~ 0.19 at intermediate lattice momenta aq around 
1, i.e. close to the scaling prediction of kz — 0.2, with a 
tendency to decrease towards a value of 0.16 at large aq. 
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FIG. 1: Gluon propagator in d = 2 for diflerent lattice sizes. 



FIG. 3: Ghost dressing function in d = 2. 
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FIG. 2: Finite-volume behavior of the gluon propagator at 
zero momentum in d = 2. Also shown are the corresponding 
bounds from (32i] . 
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FIG. 4: 'Local k\ i.e. extracted from qi and qi+s from the 2D 
gluon {k.z) and ghost (k) data, via Eqs. (|12a|) and ()12b|l . 



At this point, this would all still be quite consistent with 
the scaling solution [11,], which has also been concluded 
from previous lattice studies at finite coupling (65j . 

Our results for the ghost propagator, however, show 
much more significant deviations from this behavior. 
In particular, it appears that they cannot be recon- 
ciled with the infrared scaling hypothesis, as widely ex- 
pected to hold at least for the two-dimensional theory 
Ini [H m, 0, m [el in order to draw this con- 
clusion, it is first of all crucial to have a sound statisti- 
cal basis. In that respect it made a dramatic difference 
that we have employed a plane- wave method in obtaining 
the ghost propagator from the inversion of the Faddeev- 
Popoy operator, instead of the point-source method used 
in [61]. The two methods are compared and discussed in 
more detail in Appendix |21 

Fitting the ghost propagator, shown in Fig. [31 with 
a naive power-law ansatz at intermediate and large mo- 
menta, we extract a value of around k = 0.37. This is 
significantly larger than the expected value of 0.2 pTl.l65j| 
and than the corresponding scaling exponent from the 



two-dimensional gluon propagator in this range. More- 
over, the K from the ghost propagator continues to in- 
crease as larger and larger lattice momenta are being used 
for the fit. A good way to illustrate this, and thus the ab- 
sence of a stable scaling region for the ghost propagator 
in c? = 2 dimensions, is to extract the value of k locally, 
from just two momenta (i.e., momentum norms) via 



local 



local 



1 [ log Jgi(g.,+^) - logDgijq,) 
4 [ log Qi+s - log Qi 

log£)gh(g») - log j:)gh(g,+A-) 

log Qi+s - log Qi 



+ id-2) 



- 2 



(12a) 
; (12b) 



see Fig. H] The statistical errors have been estimated us- 
ing a bootstrap analysis.^ If anywhere, the locally defined 



^ A similar analysis was performed in three dimensions (see also 
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FIG. 5: Gluon propagator in d = 3 for various lattice sizes. 
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FIG. 6: Naive extrapolation of the gluon propagator in d = 3 
at zero momentum to infinite lattice size. The decoupling 
branch persists in this limit. The bounds from 3^ are also 
indicated. 



exponent levels around relatively small lattice momenta 
of aq around 0.1 at a value which is several standard de- 
viations below the gluon exponent there, and it continues 
to grow from there on towards larger momenta. Our re- 
sults are clearly at odds with results at finite coupling 
on moderately sized lattices [63| supposedly showing a 
scaling behavior in agreement with predictions [ill [l3 | . 



B. Three dimensions 

The gluon propagator in d = 3 is again relatively well- 
behaved. It shows a decoupling branch at small lattice 



our next subsection) already in [6lll with the important differ- 
ence, however, that the Faddeev-Popov operator was there in- 
verted on a point source. 




aq 



FIG. 7: 'Local k', as in Fig. |4] but for three dimensions. 



momenta aq and a behavior resembling a scaling branch 
at large aq, see Fig. [5] While the decoupling branch 
still shows quite significant finite-volume effects, these 
are clearly not strong enough to permit a vanishing gluon 
propagator at zero momentum in the L — >■ cxi limit, see 
Fig. [6] 

As in two and four dimensions, the branch at large aq 
is consistent with a power law with an exponent kz close 
to the value predicted by functional continuum meth- 
ods, i.e., Kz ~ 0.35 versus the predicted 0.3976 pH. 
and in agreement with [6l|. At large lattice momenta, 
the power-law exponent slightly decreases again, which 
is also evident from its local definition, see Fig. [71 

At the same time. Fig. [7] illustrates that things become 
murky again when analyzing the ghost propagator for po- 
tential power-law windows. Qualitatively, the situation 
is similar to that in d = 2 dimensions, the locally defined 
exponent continues to increase with increasing momenta 
and the ghost propagator cannot reasonably be inter- 
preted as obeying a scaling behavior in any momentum 
range. The result for the ghost propagator itself is given 
in Fig. |S] which includes a fit over a maximal data set at 
large lattice momenta (extending up to the largest ones) 
such that x^/ndf < 1. The general conclusions are much 
the same as in the d = 2 case above. 



C. The max-B gauge 

The main aspect of the results that prevents a straight- 
forward interpretation is the behavior of the ghost propa- 
gator which is most sensitive to the Gribov ambiguity. It 
therefore seems to be a natural question to ask what the 
influence is of this residual gauge freedom of the Landau 
gauge which can be used to change the ghost propaga- 
tor as explained in Sec. IIIII Here we deliberately choose 
the max-i? Landau gauge of Sec. IIIII in order to demon- 
strate the maximal effect. These are the first calculations 
at /3 = in this gauge. In the strong-coupling limit one 
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FIG. 8: Ghost dressing function for a 64^ lattice. The dashed 
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FIG. 9: Ghost dressing function in the max-_B gauge for dif- 
ferent numbers of Gribov copies taken into account. 



generates random gauge orbits for which the typical num- 
ber of Gribov copies is extremely large. Moreover, this 
number is expected to grow exponentially with the num- 
ber of lattice sites iV*. It will therefore be practically 
impossible to find even a significant fraction of all copies 
on large lattices. We therefore restrict the analysis to rel- 
atively modest lattice sizes such as 48^ and 20^. It would 
have been possible of course that the results saturated 
rather quickly with the number of Gribov copies taken 
into acco unt, as it appeared to be the case in calculations 
at P > [23,|30|. This does not happen here, however, 
at least not for the up to 600 copies per orbit analyzed 
on these lattices. As one can see below in Fig. [121 in 
particular the low-momentum values of the ghost propa- 
gator continue to grow as more and more Gribov copies 
are included in the max-B gauge. 

Our results for the ghost propagator in d = 3 dimen- 
sions are shown in Fig. [HI We observe a strong increase 
in the ghost dressing function with the number of Gri- 
bov copies ricopy included per orbit in the search for the 



FIG. 10: Effect of the max-B-gauge on the lowest nontrivial 
(i.e., non- vanishing) FPO eigenvalue as a function of included 
Gribov copies. Entries in the histogram correspond to the 
lowest eigenvalue on different configurations for the Gribov 
copy with the most divergent ghost propagator. 
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FIG. 11: Strong correlation between lowest FPO eigenvalue 
and ghost propagator at lowest non- vanishing momentum (in 
the max-B gauge). 



max-_B gauge. The observed enhancement in the ghost 
dressing function at low momenta comes along with a 
significant shift in the spectrum of the lattice Faddeev- 
Popov operator (FPO), as one would expect. This is 
illustrated on an even smaller lattice in Fig. [TU] The 
strong correlation between the lowest Faddeev-Popov 
eigenvalue Aq, omitting the three trivial zero eigenmodes, 
and the ghost propagator at lowest non-vanishing mo- 
mentum is seen in Fig. [Tl] (see also [75i] ) . 

It seems surprising that the ghost propagator in the 
max-_B gauge appears to show no sign of saturation with 
the number of Gribov copies taken into account, even on 
these small lattices and in both, two and three dimen- 
sions, as illustrated in Fig. [T^l It appears not likely that 
the ghost propagator should increase without bound as 
the number of copies is further increased in a finite vol- 
ume, but we are obviously still far from having sampled 
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FIG. 12: Ghost dressing function vs. inverse number of Gribov copies for the max-B gauge. Left: for d — 2 dimensions and a 
48'^ lattice. Right: for d = 3 dimensions and a 20^ lattice. Lines are drawn to guide the eye. 
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FIG. 13: Local kz resp. k from the max-B gauge. 



sufSciently many to conclude that we could reach a stable 
limit. We conclude that Gribov copies do have an enor- 
mous effect on the low-momentum behavior of the ghost 
propagator. This is also reflected in the corresponding 
k's as shown in Figure [T^ 

It should be noted that by choosing instead of the max- 
imum value of the ghost propagator a finite value for B, 
it is possible to generate any ghost propagator in be- 
tween the one of minimal Landau gauge and the one in 
the max-i? gauge. In fact, it is also possible to select a 
ghost propagator which is smaller than the one in mini- 
mal Landau gauge 0^ . 

However, our primary interest in employing the max-i? 
gauge is to investigate a potential fixed-point behavior in 
the running coupling, an issue to which we now turn. 



1 

aq 

FIG. 14: 'Local k\ as in Fig. |3]but for d — 4 dimensions; see 
[43 | for the corresponding propagator results. For the ghost 
propagator, not all momenta eligible under the cylinder cut 
have been used here. The exponent nz from the gluon data 
is virtually constant except for the smallest lattice momenta. 
Error bars are relatively small and not shown. 



D. Running coupling 

Let us emphasize that an increase of k from the ghost 
propagator towards large lattice momenta is actually also 
apparent from the d — 4 data presented already in [42I - 
1441 ■ as shown explicitly in Fig. [TH But it is important 
to notice that in this case, unlike in d = 2 and d = 3, 
Fig. [T21 K from the ghost data does not rise above kz, 
i.e., the ghost exponent extracted from the gluon data. 
This manifests itself in a different behavior of the effective 
running couphng As reported in m,|43|, in d = 4, as 
grows monotonically as a function of aq and approaches 
a value ofc ~ 4 from below, which is close to the value 
Q!c ~ 4.46 predicted in p^ . 

The behavior we observe here in minimal Landau 
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FIG. 15: ae in d = 2, L = 288, and d = 3, L = 64, as well as 
from various lattice sizes in d = 4. 



gauge in d = 2 and d = 3 differs surprisingly: The be- 
havior at large lattice momenta is notably different from 
the four-dimensional case, see Fig. [151 This is obviously 
a joint consequence of the behavior of k and kz (Figs. 0] 
respectively [7]) and the fact that from Eq. © 

aeOc(<72)2(-.-). (13) 

In c? = 2 and 3, where the local k grows monotonically 
and intersects kz at some aq, see Figs.|3]and[71 ae reaches 
a maximum near the respective lattice momentum, as 
implied by Eq. (I13p . Comparing the results to the pre- 
dictions from : 12ij , we find the maximum of is approx. 
85% of the predicted (constant) value of etc in d = 2 and 
« 70% in d = 3, instead of « 90% in d = 4. This refers 
to the prediction for ac given the standard values of k, 
i.e. 0.3976 ... in d = 3 and 0.2 in d = 2. The situation 
slightly improves if the value of the predicted ac is calcu- 
lated instead with the k values actually observed on the 
lattice. In d = 3, this amounts to choosing k = 0.35, see 

Fig. nil 

While the behavior of Oe in d = 2 and 3 is surprising 
and not amenable to an immediate explanation, it is in- 
teresting to note that the maximum of as{aq) shows a 
monotonic behavior from d = 2 to d = 4, i.e., it is shifted 
to larger aq. This follows the general pattern observed 
also at finite /?: the window in which a scaling-like behav- 
ior is observed extends to much smaller momenta in two 
dimensions than in three dimensions, before the effective 
masses begin to dominate. In four dimensions, finally, the 
window becomes so small that almost no scaling-like be- 
havior can yet be observed. The same mechanism seems 
here to shift the observed maximum to ever smaller val- 
ues of aq. 

A possible explanation for these results would be that 
at large aq discretization effects are visible while at small 
aq a Gribov copy effect causes the deviation of ae from 
the continuum solution {a,, — ac — const). 

The second possibility motivates us to investigate the 
effect of the max-_B gauge with a growing number of 
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FIG. 16: max-B gauge: Very strong change in the running 
coupling, especially at small lattice momenta, as the number 
of Gribov copies increases. Top: for d = 2 dimensions and a 
48^ lattice. Bottom: for d = 3 dimensions and a 20^ lattice. 



copies on ae. We have used up to 600 Gribov copies per 
configuration on a 20'^ lattice. The effect on the ghost 
propagator is strong, as is apparent from Fig. [121 (right) 
in the previous subsection. We observe no 'saturation'. 

By virtue of Eq. ([9]), this effect induces an even 
stronger change in the running coupling (see Fig. I16p . 
which is far from being balanced by a similarly strong 
change in the gluon propagator. In fact, it was already 
stated in that the effect of choosing the max-B gauge 
on the gluon propagator is relatively small for the em- 
ployed lattice sizes. 

Since the resuh overshoots the prediction this implies 
that it is possible to select a value of B such that actually 
an infrared finite coupling could be generated. Such a 
choice leads to essentially the same behavior, as is seen 
when artificially the number of copies is restricted to ~ 70 
in d = 3. This is equivalent to the statement that the 
locally extracted values of k and kz are consistent within 
errors in this region, see Fig. [T3l 

There is currently no fully satisfactory explanation for 
the behavior of the running coupling and the observed 
over-scaling. Even though a ricopy — ^ oo extrapolation 
is not at all feasible from the current data, the result 
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in this limit on a 20^ lattice will be at odds with the 
continuum prediction for ac- One possible solution might 
be to consider larger lattices. 

We have done similar simulations for the two- 
dimensional case, again initially on a small lattice 
(48^). The qualitative picture is similar to the three- 
dimensional case, see Fig. as a scaling region is also 
absent. As noted above, this is surprising. In particular, 
when considering the scaling behavior observed in two 
dimensions [s^, 0, 113, HE] at volumes where only very 
few Gribov copies are present [i^ . 

Regarding the presence of possible discretization ar- 
tifacts at large momenta, one possible concern in this 
context is the lattice breaking of rotational invariance at 
/3 — 0. The fact that the dressing functions on the lattice 
are invariant under the lattice isometrics does not imply 
that they are also invariant under the Euclidean 0{d) 
symmetry in d dimensions. This implies that the dressing 
functions depend on the orbits generated by the lattice 
isometrics, but there are different such orbits that corre- 
spond to the same 0{d) orbit classified by the magnitude 
of momentum (see, e.g., [76]). These differences are of 
higher order in the lattice spacing and would vanish in 
the continuum limit. There is no suppression of these dif- 
ferences in the strong-coupling limit, however, and they 
will therefore show up in differences of dressing functions 
evaluated along different momentum directions. We as- 
sess these effects of the lattice breaking of rotational in- 
variance in Appendix |B] They are not as dramatic as 
the Gribov effects, but the tendency of the fixed-point 
couplings to decrease at large lattice momenta is within 
the systematic uncertainty due to these effects. 



V. CONCLUSIONS 

We close with a summary of our main results and a dis- 
cussion of the consequences. In the present work we have 
qualitatively extended previous studies on the strong- 
coupling limit /3 = of Landau gauge SU(2) Yang-Mills 
theory on the lattice. In two and three dimensions we 
have shown that the Gribov ambiguity is rather strong. 
It especially affects the ghost propagator which neither 
uniquely shows a scaling nor a decoupling behavior for 
small lattice momenta. Indeed, subject to the number 
of gauge copies used, it even shows over-scaling for small 
and intermediate momenta. Our /3 = results also in- 
dicate that over-scaling can be achieved for any finite 
lattice if enough gauge copies are taken into account. 
However, in the continuum over-scaling is prohibited by 
the uniqueness proof put forward in [IEl [igI . [70I |. It is 
not part of the one-parameter family of solutions 0, H^l 
and might hint at the scaling solution as the limit of the 
over-scaling lattice Landau gauge fixing. 

In any case, the above results are in marked contrast 
to the results in previous works where only a deco upl i ng- 
type solution was obtained at low momenta, see |42h441 
|61 |. Evidently the interpretation of the /? = results 



needs to be seriously reconsidered. 

At large lattice momenta, the ghost propagator does 
not naively exhibit scaling, as its local scaling coefficient, 
K, rises monotonically. This is neither a scaling behavior, 
nor is it a decoupling-type behavior. Indeed, a logarith- 
mic fit appears to be possible, and has been related to 
similar fits for the decoupling solution in the continuum 
(61i] . Note, however, that this similarity is superficial be- 
cause the prefactor of the logarithmic term has the wrong 
sign. A naive extrapolation to finite /3 and the infinite- 
volume limit would lead to a pole at some finite momen- 
tum. On the other hand, the gluon propagator shows 
approximate scaling and hence a qualitatively different 
behavior from the ghost. Such a different behavior would 
not be expected for a decoupling-type solution. A pos- 
sible explanation for these surprising results are strong 
but not unexpected discretization effects associated with 
the unphysical /? — limit at large aq, see Appendix IB] 

Our results are best summarized in Fig. ll6l for the run- 
ning coupling. It is the only RG-invariant product of 
the two dressing functions and is rather sensitive to the 
global aspects of the gauge fixing. Its behavior at small 
lattice momenta confirms that the Gribov ambiguity is 
much stronger than expected on the basis of previous re- 
sults. Extrapolating to finite /3 in four dimensions, we 
conclude that a resolution of the Gribov ambiguity re- 
quires far more computational power than presently em- 
ployed 32-34, 36, 37, 77]. In particular this implies that 
the present data cannot be used to exclude any subset of 
infrared solutions. It is all within the uncertainties due 
to the Gribov ambiguity. Indeed, we are confident that 
our data together with other lattice and continuum re- 
sults 0,[Il|2l,[2l,|4l|43,[53-[5i[2l[7i can be unified in 
an overall picture of Landau gauge Yang-Mills theory as 
put forward in 0, 113, [s^] • Its further understanding 
would certainly provide more insight into the confining 
physics. 
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FIG. 17: Ghost dressing function in d = 2. In order to draw 
conclusions about the behavior at large lattice momenta, the 
plane-wave source method is clearly preferable to the point- 
source method: only 76 meas. on 288^ (plane -wave source) vs. 
~ 37000 on 280^ (point source) yield much less fluctuations. 



FIG. 19: The ghost dressing function in three dimensions us- 
ing the point-source method. The numbers of independent 
configurations are given in the legend. For an explanation of 
the wiggles at large momenta see [t^ . 
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FIG. 18: 'Local At' for a point source from 37000 measure- 
ments and calculated from (cylinder cut) momentum norms 
that are 5 = 20 indices apart. Compare to the result from a 
plane-wave source with 76 meas. at = 288^ in Fig. |4l 



Appendix A: Plane-wave vs. point source 

We expand here the study of [73i], which showed that 
the method to obtain the inverse of the Faddeev-Popov 
operator to calculate the ghost propagator has quite an 
impact on the statistical accuracy of the calculation. Es- 
pecially, for extracting information about the local be- 
havior it is advantageous to employ a plane- wave source 
[40| instead of a point source [t^ for the conjugate gradi- 
ent method used to invert the operator. This means that 
the Faddeev-Popov operator is inverted on a vector of 
plane waves, s^'^(A:) = S°'^e'^''''" , instead of a point source, 
Sfc" = ^"'(^-,0 - ^/V), with k-x = 2TTj:^ k^xJL^. 

The contrast is illustrated in Fig. [T71 which already 






20', n =550 

copy 


A 


n =300 

copy 




n =80 

copy 


□ 


n =20 

copy 


o 


"copy"' 


X 


80^ n =20 

copy 


+ 


n =1 

copy 





FIG. 20: Finite-volume effect for minimal Landau gauge and 
for max-B gauge with various numbers of copies. Point-source 
method. 



suggests that the point source makes a power-law fit at 
large lattice momenta, and thereby the identification of 
a scaling region with an unique exponent k, consider- 
ably more difficult. This statement can be made more 
precise, e.g., by performing power-law fits at large lat- 
tice momenta and comparing the respective values of 
X^/ndf. The reason is that a plane-wave source aver- 
ages over all possible inversion points, and therefore pro- 
vides a volume-factor less noise than a point source. The 
drawback of a plane-wave source is that an inversion is 
required for every momentum value. 

This problem is further illustrated in Fig. [18] for a lo- 
cally defined k, in contrast to Fig. S] (plane- wave source). 
Of course, in the limit of infinite statistics, both meth- 
ods will yield the same result. Using a point source to 
determine the value of k at large momenta, a value of 
0.31 would be obtained. 
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FIG. 21: Effective running coupling in d = 2 for two different 
choices of lattice momenta. 



Similar considerations apply in all dimensions. The 
result for the point-source method in three dimensions is 
shown in Fig. [THl The fluctuations are greatly enhanced 
compared to the plane- wave result shown in Fig. [5] Even 
with N € C(IO^) configurations, they are still larger than 
in the latter case. Correspondingly, the region where a 
power-law fit at large aq with x^/ndf < 1 is possible is 
much smaller. 

There is also a note in conjunction with the max-i? 
gauge employed. The ghost propagator using the point 
source exhibits much larger statistical fluctuations. Con- 
sequently, when triggering on large fluctuations, as when 
implementing the max-i? gauge, it is likely that in a 
sample of finite statistics the ghost propagator will be 
larger in the max-i? gauge than when using a plane-wave 
source. Due to its inherent averaging over the lattice 
the plane-wave source method smooths these exceptional 
fluctuations in many more cases. Only when volume- 
times more statistics have been used for the point-source 
method than for the plane-wave source method, a similar 
smoothing effect can be expected to take place. Again, 
in the limit of infinite statistics both methods will then 
yield the same result. A measure of this, aside from the 
ghost propagator directly, is the comparison of average 
vs. median. At the present statistics, in both cases both 
measures significantly deviate but the difference is much 
larger for the point-source method than for the plane- 
wave source method. Also, when using the median in- 
stead of the average the results for both methods are 
more similar. 

As an example, the results at finite statistics for the 



point-source method are shown in Fig. 1201 Comparing to 
Fig. E] from the main text displays exactly this point. 



Appendix B: Discretization effects 

As noted above, we have usually applied a cylinder cut 
to the momenta. In order to obtain a rough estimate of 
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FIG. 22; Effective running coupling in d = 3 for tliree different 
choices of lattice momenta. 



possible discretization artifacts, we have also calculated 
the running coupling for the opposite case of on-axis mo- 
menta, ki = kSij for all fixed j. Interestingly, the effec- 
tive running coupling thus obtained is compatible with 
a constant value at large lattice momenta close to the 
continuum prediction Tz| , see Fig. [5T] These data have 
been produced with the point-source method, which suf- 
fices to clearly observe the qualitative effect. Data at 
diagonal momenta, which form a subset of the momenta 
surviving the cylinder cut, are shown for comparison. In 
d — 3, the difference between the choices of momenta is 
again clearly visible, see Fig. [511 though less pronounced 
than in the lower-dimensional case. This suggests an even 
weaker effect for d = 4, where indeed the running cou- 
pling has been found to tend towards a constant value 
at large aq [H, These early results underline the 

serious discretization problem faced at large aq at (3 = 
and call for further investigations along these lines, e.g. 
in max-i? gauge. 
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